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Abstract 

A survey is presented of several extrema principles of 
energy dissipation as applied to problems in fluid mechanics. An 
exact equation is derived for the dissipation function of a homo- 
geneous, isotropic, Newtonian fluid, with terms associated with 
irreversible compression or expansion, wave radiation, and the 
square of the vorticity. By using entropy extrema principles, 
simple flows such as the incompressible channel flow and the 
cylindrical vortex are identified as minimal dissipative distribu- 
tions. The principal notions of stability of parallel shear flows 
appear to be associated with a maximum dissipation condition. 
These different conditions are consistent with Prigogine’s clas- 
sification of thermodynamic states into categories of equilib- 
rium, linear non-equilibrium, and non-linear non-equilibrium 
thermodynamics; vortices and acoustic waves appear as exam- 
ples of dissipative structures. The measurements of a typical 
periodic shear flow, the rectangular wall jet, show that direct 
measurements of the dissipative terms are possible. 
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speed of sound in undisturbed fluid 

A,B,C 

integration constants 

djk 

deformation rate tensor 

h 

parallel channel height 

Ii. h 

dissipation integrals (also I T , 1^) 

k 

coefficient of heat conduction 

K 

vortex strength 

L 

length of vortex 

Mj 

Mach number in the j direction = uj/ao 

P 

pressure 

r 

radial distance 

r o 

vortex core radius 

t 

time 

u 

stream wise velocity component 

u, Uj 

velocity vectors 

U0 

tangential velocity component 
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U 0 maximum nozzle exit velocity 

W parallel channel width 

x stream wise direction 

x, Xj position vector 

y direction normal to parallel channel wall 

5j k Kroneker delta 

O viscous dissipation function 

incompressible dissipation function 
<I> R radiative component of dissipation function 

X second coefficient of viscosity 

\i viscosity 

v kinematic viscosity 

p density 

co, C 0 j k vorticity 

( ) time average of quantity in parentheses 

( )* time dependent fluctuation of quantity 

in parentheses 


L Introduction 

The successful application of the "principle of least 
action" to problems in the fields of particle dynamics and quan- 
tum mechanics has inspired several efforts to extend the method 
to other areas, including thermodynamics and hydrodynamics. 
The relevance of variational principles to these areas has been 
studied relatively recently, and hence there is significantly less 
acceptance of any general framework of extrema principles to 
them. Several writers [1,2] conclude that variational principles 
are not applicable to the field of thermodynamics. Onsager 
examined coupled irreversible processes near equilibrium and 
derived a statement of minimum entropy production for these 
conditions [3], Ziegler [4] proposes that rates of entropy pro- 
duction and energy dissipation arc maximized for processes far 
from equilibrium, while Biot [5] supports Onsager’s conjecture 
of minimum dissipation rate for near-equilibrium processes. 
Tadmor [6] has also defined a minimum entropy principle in the 
context of stable numerical solutions of the gas dynamics 
equations. 

In the field of hydrodynamics, Helmholtz [7] proved that 
flow with negligible inertia is characterized by a minimum of 
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viscous dissipation. Herivel and Lin [1] have derived iirota- 
tional and rotational forms respectively of the conservation 
equations from a variational principle, as related by Yourgrau 
and Mandelstam [1], For thermal convection and turbulent 
flow in channels, Malkus [8,9] proposed a theory based on 
arguments that the smallest turbulent scales are marginally sta- 
ble, and that the mean and turbulent characteristics of the flow 
are determined by the condition that the energy dissipation is 
maximized. Objections to Malkus’s theory have been raised in 
several areas, among them that extrema of the dissipation func- 
tion are not generally compatible with conservation laws, and of 
vagueness in the application of both extremizing conditions and 
stability criteria. Reynolds and Tiederman [10] found poor 
agreement between a eddy viscosity simulation of turbulent 
channel flow and Malkus's prediction of the association of 
maximum dissipation rate with neutral stability of the mean 
flow. Nihoul [11] reexamined Malkus‘s theory with regard to 
Liapounov stability and described an associated principle of 
minimum Reynolds number. The problem of convective heat 
transfer, also considered by Paltridge [12] and Busse [13], can 
be characterized by maximum heat transfer in some 
circumstances. 

Prigogine [14] has classified entropy production theo- 
rems according to a generalized sequence of stable thermody- 
namic states. Thermodynamic equilibrium, characterized by the 
absence of gradients of state or kinematic variables, is in a state 
of maximum entropy and zero entropy production. Linear 
nonequilibrium processes, studied by Onsager [3], are associ- 
ated with minimum entropy production. The linear description 
applies if the generalized thermodynamic fluxes (such as heat 
flow or fluid deformation) are linearly related to the generalized 
thermodynamic forces (such as temperature gradient or viscous 
stress). Many near-equilibrium processes are accurately charac- 
terized as linear. Entropy production is not necessarily mini- 
mized for nonequilibrium processes which arc nonlinear or far 
from equilibrium. Nonlinearity is generally associated with the 
generation of chaotic states from initially determinate conditions, 
such as the transition of laminar flow to turbulence. Prigogine 
[15,16] has demonstrated that nonlinearity and processes which 
are far from equilibrium may give rise to global organization 
from an initially random field, and associates such processes 
with the evolution of dissipative structures. Prigogine describes 
several general characteristics of nonlinear systems, such as limit 
cycles and period-doubling Feigenbaum sequences, which are 
familiar aspects of high speed and transitional flows. Lugt [17] 
has commented on the evolution of discrete vortices in shear 
flows as an example of these dissipative structures. 

II. Dissipation functions for a Newtonian fluid 
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Hinze [19] and Tennekes and Lumley [20] obtain similar 
decompositions of the time-averaged dissipation function, which 
are valid for homogeneous turbulence. The last term in (3) 
appears in several aeroacoustic wave equations, such as that 
derived by Phillips [21]. Home and Karamcheti [22] have 
derived a different formulation of the dissipation function by 
combining the terms of the equation: 

V- + = 0 (4) 


where: 


3p 

at 


3p 


3u i 


3u. du. -s n ax' 

F = p — !+pu i + 1 

H at H j ax. ax. ax t 

i * j 


(5) 

( 6 ) 




au. a u 


ax. 


ax, 


au 

+ ValT 


M = F = 0 for a fluid free of sources of mass or momentum. By 
expanding the terms of equation (4), we obtain: 
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The present discussion examines the dissipation function 
for a homogeneous Newtonian fluid, which Ziegler [4] writes 
as: 
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Jeffreys [18] rewrites (1) as: 
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By cancelling terms as shown, the remaining terms in equation 
(4) appear as follows: 


rate. For parallel channel flow, the mass and momentum equa- 
tions reduce to: 
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By substituting equation (7) into equation (3), the 
desired form is obtained: 
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This exact expression accounts for three important 
mechanisms of energy dissipation in real flows: (1) irreversible 
expansion or compression, (2) generation and radiation of sound 
or shock waves, and (3) generation of vorticity. Since the mass 
and momentum conservation equations have been used in the 
derivation of (8), extrema of this function will satisfy those 
equations. Formally, functions of the conservation equations 
have been adjoined to the conventional expression for the dissi- 
pation function similar to the manner of imposing functional 
constraints in the solution of problems of geodesics. The radia- 
tive terms can be further manipulated to yield a convective wave 
operator acting on the density field: 
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For this case, the term in the dissipation function associ- 
ated with the Laplacian of the pressure is zero, as is seen by 
taking the x-derivative of the momentum equation and by using 
the continuity equation: 
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The volume integral of the dissipation function is now 
given by: 

> I dv=p|(o 2 dv = pxw £ dy 0°) 


where x and w are the length and width, respectively, of the 
channel. The flow rate through the channel is given by: 


V c = Judy 

Possible extrema of the integral (10), subject to constant 
channel flow rate, arc found as solutions of the Eule. Lagrange 
equation: 
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where is the Mach number in the k direction. The first term 
occurs in Lighthill’s [24] acoustic stress tensor, while the con- 
vective wave operator for the density is identical to the convec- 
tive wave operator for acoustical disturbances in shear flow. 
The presence of non-dissipative outward radiating acoustic 
waves from an aerodynamic source region requires dissipation 
of energy from the source region. 


Here, a is a Lagrange multiplier adjoining the integrand of the 
flow rate constraint to the extremizing function. For this exam- 
ple, evaluation of the Euler-Lagrange equation yields the 
expression: 


For incompressible, uniform-density flows, the expres- 
sion (8) for the dissipation function reduces to: 


*, = -—? P+P® 


This expression may be used to examine low-speed flows for 
velocity distributions which result in extrema of the dissipation 


Hence, an extremum of the dissipation is found for the linear- 
parabolic profile: 

u(y) = Ay 2 + By + C 

This profile corresponds to a minimum of the dissipation func- 
tion, according to Legendre's test, since: 
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Both solutions represent minimum dissipation conditions, as in 
(13) the channel flow. Integrated contributions to the overall 
dissipation rate are given in the following table: 


We seek extrema of the integral: 

| <t>dv = - jiL J V 2 (p) rdr d9 + ji2kL £ a> 2 rdr 

where L is the length of the vortex. The first integral on the 
right side vanishes, since from equation (1 1): 

J 7 a (p)rdr= f }£(u 2 Jrdr=u 2 | 

00 0 

and 

u e (°) = u # («) =0 
for any realizable flow. Then: 

J*dv = 2 *L J (^i) + £ rdr (14) 

0 

Extrema of this integral are obtained as solutions of the Euler- 
Lagrange equation: 
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As can be seen, the net dissipation in the rotational core is zero, 
while the irrotational outer flow has a net positive dissipation 
due to the Laplacian of the pressure. The tangential acceleration 
is zero throughout the flow, except at the boundary at r Q . A 
steady tangential stress is required at this location (as provided 
by a thin rotating cylinder, for example) to maintain a steady 
flow. The net dissipation in the outer region is identical to the 
torque power absorbed by the rotating cylinder, since the shear 
stress at the cylinder wall is given by: 



The two example considered are steady flows, since the 
Euler-Lagrange equations did not include time dependent varia- 
tions. The two examples also involve no convective accelera- 
tions, and therefore do not demonstrate an extension of the 
minimum dissipation principle beyond Helmholtz's original 
statement for flows with negligible inertial forces. For these 
flows, the minimum dissipation condition appears to be con- 
sistent with a linear, near-equilibrium process. It should be 
possible to extend the method to boundary layers and asym- 
metric flows where convective accelerations are important 
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IV, Dissipation considerations in unst eady, two-dimensional 
flows 

The application of extremum conditions of entropy pro- 
duction or dissipation rate to unsteady flows may be facilitated 
by considering the mean and fluctuating components of the dis- 
sipation function: 

p = P+p'; co = S3 + co' 

Oj =- -^£v 2 (p + p') +p(E + coO 

+ n[(®) +(<o') ] (18 ) 
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The Laplacian of the pressure was computed from 
the divergence of the incompressible, viscous, 2-D momentum 
equation: 


We may use this expression to interpret the small and 
large disturbance motions in two-dimensional flow. A large 
number of parallel shear flows such as the boundary-layer, the 
free-shear-layer, and the wall-jet have been successfully ana- 
lyzed for stability to small disturbances via the Orr-Sommerfield 
equation. The stability characteristics for these flows are found 
to be strongly associated with the mean velocity profile of the 
flow. For velocities which exceed a critical Reynolds number, 
small disturbances of the frequency corresponding to the maxi- 
mum amplification rate are predicted to grow exponentially with 
downstream distance until nonlinear effects limit the growth. 
Within the small -disturbance region, the mean velocity profile 
remains unchanged, and the small disturbances take the form of 
convecting vortical motions. Experimental studies of various 
flows confirm that the frequencies predicted for maximum- 
vortical-disturbance growth rate correspond to the observed 
frequencies of unforced fluctuations. From equation (18), it 
appears that the maximum vortical-disturbance growth rate con- 
dition corresponds to a maximum dissipation condition, if the 
contribution from the pressure term is negligible and the contri- 
bution from the mean vorticity term is independent of distur- 
bance frequency. 

Figure 1 illustrates the periodic vortex motion of a two- 
dimensional wall jet flow. The wall jet develops from an 
initially parabolic profile rather than a self-similar profile, as is 
the case in the parallel shear flow. Woolley and Karamcheti [25] 
have shown that the stability characteristics of nonparallel shear 
flows are closely related to those of parallel shear flows. The 
short wall was associated with the generation of tones by the jet, 
which significantly reduced the small disturbance region by 
effectively forcing the initial jet region, however the observed 
tone frequencies were approximately the same as with the long- 
wall, silent jet. The conditions for this jet are: U m = 13.85 
m/s, L/h = 7.5, where U m = the maximum velocity of the 
parabolic exit profile, L = wall length, and h = nozzle width. 
Figure 1 (a) shows a phase- averaged Schlieren visualization 
record obtained by lightly heating the subsonic nozzle flow. 
Figure 1 (b) shows the phase-averaged velocity field, referenced 
to the convecting vortices, and Fig. 1(c) depicts the 
corresponding vorticity field obtained with a central-differenced 
curl of the velocity field. The vorticity is normalized with 
respect to the maximum exit velocity and the nozzle width. The 
measurements were obtained from a single x-wire velocity probe 
by sampling the probe output at regular phase intervals as 
determined by a fixed pressure-transducer in the wall [23]. 
These measurements were further processed to obtain estimates 
of the mean and fluctuating components of the dissipation field. 

Relative contributions to the dissipative structure of the 
wall jet were estimated by integrating the field variables in the 
cross- stream direction. We define: 





In order to estimate the measurement error, the mean square of 
the divergence of the velocity was computed: 



Figure 2 shows the distributions of Ii , h , I 3 , Idiv . 
and It for the jet conditions described previously. The mea- 
surement error indicated by Idiv is negligible only upstream of 
the wall trailing edge, where the pressure Laplacian term, Ii , is 
also negligible. The square of the mean vorticity, I 2 , is nearly 
constant in the small disturbance region: 0 < x/h , 2 . 5 , beyond 
which it gradually decreases. The mean square of fluctuating 
vorticity, I 3 , is negligible in the initial region, then steadily 
increases in the large disturbance wall region: 2.5 < x/h < 7.5. 
Relative comparisons of the magnitudes of the pressure and 
vorticity fluctuation terms in the small disturbance region are not 
possible with these measurements, but could be made from a 
stability analysis of the flow. 

These measurements demonstrate a potential method of 
direedy measuring dissipative terms in an unsteady flow. A 
direct extension of this method may be applied to experimentally 
determine a relationship between overall dissipation and variable 
parameters, such as forcing frequency, and to exp imentally 
search for dissipation extrema. 

^■Conclusions 

An exact equation has been derived for the dissipation 
function of a homogeneous, isotropic, Newtonian fluid, with 
terms associated with irreversible compression or expansion, 
wave radiation, and the square of the vorticity. Simple flows 
such as the incompressible channel flow and the cylindrical 
vortex are identified as minimal dissipative distributions. The 
principal notions of stability of parallel shear flows appear to be 
consistent with maximum dissipation conditions on the growth 
of vortical disturbances. These observations are consistent with 
Prigogine's [14,15] distinction between stable thermodynamic 
states or processes which are cither linear and near equilibrium, 
or nonlinear and far from equilibrium. In this context, vortices 
and acoustic wave sources appear to be examples of dissipative 
structures. The measurements of a typical periodic shear flow, 
the rectangular wall jet, show that direct measurements of the 
dissipative terms are possible. 

Further experimental and theoretical analyses are 
required to determine the global validity of extrema principles 
with regard to complex flows and to determine the applicability 
of time-averaged analyses. 
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b. 



Fig. 1 Phase averaged velocity field and flow visualization, 
a) Phase-averaged flow visualization; b) velocity vectors 
(relative to convecting vortices; and c) vorticity contours. 
Conditions: nozzle width, h = 0.508 cm; nozzle aspect 

ratio = 20; wall length, L = 3.81 cm; parabolic velocity 
profile at nozzle exit with maximum exit velocity, 
U 0 = 13.85 m/sec; tone frequency, f = 600 Hz. 
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Fig. 2 Distribution of integrated dissipation components. 
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